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14  ABSTRACT 

We  have  established  Dislocation  Dynamics  (DD)  simulations  as  a  practical  method  to  compute  the  acoustic  nonlinearlity  parameter,  beta,  which 
forms  the  bases  of  the  emerging  metal  fatigue  detection  technique  of  nonlinear  ultrasonics.  Our  study  have  uncovered  mistakes  in  two  well-known 
theoretical  models  in  this  field  and  leads  to  a  deeper  understanding  of  the  dislocation  mechanism  to  acoustic  nonlinearity.  Fora  single  dislocation 
bowing  in  its  glide  plane  between  two  pinning  points  (i.e.  the  monopole  model),  our  simulations  and  analytic  derivations  show  a  strong  dependence 
of  beta  on  the  dislocation  orientation,  which  is  missed  by  the  previous  model  (Hikata  et  al.  1965).  For  parallel  dislocations  forming  a  multipole 
configuration  (as  a  model  for  the  vein  structure  in  fatigued  metal),  our  simulations  show  a  pronounced  dependence  of  beta  on  the  applied  stress 
(with  betaK)  at  zero  stress),  which  contradicts  the  previous  model  that  predicts  a  constant  beta  independent  of  stress  (Cantrell  et  al.  2001 ).  Our 
analytic  derivations  have  pinpointed  the  mistake  in  the  previous  model 
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1.  Single  Dislocation  Line  Interaction  with  Ultrasound 

We  have  studied  the  process  of  dislocation  interaction  with  an  oscillating  stress  field  in 
great  detail,  both  using  Dislocation  Dynamics  (DD)  simulations  and  analytic  models.  A 
classical  problem  is  illustrated  in  Fig.  1,  where  the  bow'ing  out  of  the  dislocation  line 
(pinned  at  both  ends)  produces  a  strain  that  is  non-linear  with  the  applied  stress.  This 
gives  rise  to  acoustic  non-linearity  /?,  which  can  be  measured  experimentally  from  the 
amplitude  of  second  harmonics. 


Figure  1.  Schematics  of  a  dislocation 
segment  with  ends  pinned  by  obstacles 
vibrating  under  the  stress  field  of  an 
ultrasound. 


Previous  studies  have  assumed  a  simple  model  for  dislocations  —  the  dislocation  is 
represented  as  a  string  with  a  constant  line  tension  (Hikata  et  al.  J.  Appl.  Phys.  36,  229, 
1965).  The  model  predicts  linear  dependence  of  beta  on  the  applied  stress  a  (for  small  ct) 
in  the  case  of  a  single  dislocation  pinned  at  two  end  points.  Our  explicit  Dislocation 
Dynamics  simulations  showed  that  the  acoustic  non-linearity  parameter  beta  depends 
critically  on  the  orientation  of  the  dislocation  and  material’s  Poisson  ratio  (see  Fig.  2). 
The  edge  and  screw  orientation  becomes  equivalent  only  in  the  limit  of  Poisson's  ratio 
going  to  0.  But  for  most  materials,  Poisson’s  ration  is  around  0.3,  in  which  case  the 
behavior  of  edge  dislocations  is  very  different  from  that  predicted  by  the  simple  line 
tension  model.  At  the  same  time,  the  microstructure  of  fatigued  metals  is  dominated  by 
edge  dislocations.  This  means  the  previous  understandings  of  dislocation  microstructurc 
-  ultrasound  interaction  needs  to  be  drastically  altered. 

We  found  that  the  deviation  from  the  explicit  DD  simulations  and  constant  line-tension 
model  is  mainly  due  to  the  dependence  of  dislocation  line  tension  on  line  orientation.  By 
allowing  the  line  tension  to  depend  on  line  orientation,  we  have  obtained  analytic 
expressions  for  the  acoustic  non-linearity  that  agrees  well  with  our  DD  simulations. 
These  results  are  described  in  a  manuscript  currently  under  review  in  Acta  Matcrialia  [7], 
This  work  corrects  an  important  error  in  previous  predictions  of  acoustic  non-linearity 
and  this  the  first  time  fi  is  calculated  explicitly  from  DD  simulations.  For  the  first  time, 


we  show  dislocation  contribution  to  P  can  be  negative.  It  means  sometimes  dislocation 
mechanisms  can  reduce  the  measured  value  of  /}.  This  naturally  explains  the  decrease  of 
P  in  some  of  the  experimental  data  reported  in  the  original  paper  of  Hikata  et  al.  (1965), 
see  Fig.2(c),  which  was  difficult  to  explain  with  the  previous  model  that  ignored 
orientation  effects. 


Figure  2.  (a)  Predicted  stress-strain  relationship  for 
edge  and  screw  dislocations  in  an  elastic  medium 
with  different  Poisson’s  ratio  v  when  orientation 
effects  are  included,  (b)  Predicted  acoustic  non¬ 
linearity  parameter  P  based  on  the  stress-strain 
relationship  in  (a).  Notice  P  becomes  negative  in 
some  stress  regime  for  edge  dislocations  in 
materials  with  v  =  0.3.  (c)  Experimental  data  for  p 
as  a  function  of  applied  stress  in  the  original  paper 
by  Hikata  et  al.  (1965)  for  aluminum  samples  pre- 
deformed  to  different  amount  of  strains.  Notice  the 
initial  decrease  of  p. 


2.  Dislocation  Dipole  and  Multipole  Interaction  with  Ultrasound 
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Figure  4.  (a)  Acoustic  nonlinearity  p  of  a  dislocation  dipole  (left)  as  a  function 
of  applied  stress  predicted  by  DD  simulations  (stars)  and  our  analytic  solution 
(solid  line).  The  prediction  from  the  previous  model  (Cantrell  2001)  is  plotted 
as  dashed  line,  (b)  Acoustic  nonlinearity  P  as  a  function  of  applied  stress  of  a 
group  of  dislocations  forming  a  Taylor  lattice  (left)  by  DD  simulations. 


Figure  5.  Slope  of  the  curves  in  Fig.4  (i.e.  dp/dcr)  at  cr=0  as  a  function  of  the 
size  n  of  Taylor  lattice  (circles),  n  is  the  number  of  dislocations  in  the  Taylor 
lattice.  The  straight  lines  are  drawn  to  indicate  that  dp/da  scales  as  n  . 


In  the  case  of  two  parallel  dislocations  forming  a  dipole,  an  earlier  model  (Cantrell  et  al. 
Int.  J.  Fat.  23,  S487,  2001)  predicts  that  p  is  a  (positive)  constant  that  is  independent  of 
stress  a  (see  Fig.4(a),  dashed  line).  This  prediction  has  been  used  extensively  in 
subsequent  publications  to  estimate  the  value  of  P  for  configurations  containing  multiple 
dislocations.  However,  our  DD  simulations  show  that  P  has  a  pronounced  dependence  on 
stress.  Furthermore,  we  have  performed  analytic  derivations  to  pinpoint  the  mistake  in 
the  previous  solution  ---  the  series  expansion  has  not  been  carried  out  to  sufficient  order. 

For  parallel  dislocations  forming  a  multipolc  configuration  such  as  a  Taylor  lattice 
(which  is  a  model  for  the  vein  structure  in  fatigued  metal),  the  earlier  model  would 
estimate  the  value  of  P  by  multiplying  the  value  of  the  zero-stress  value  of  P  for  a  dipole 
by  the  number  of  equivalent  dislocation  dipoles  in  the  Taylor  lattice.  This  would  again 
predict  a  P  value  that  is  independent  of  the  applied  stress.  However,  our  DD  simulations 
show  that  the  p  value  at  zero  stress  is  actually  zero,  and  that  p  has  a  strong  stress- 
dependence  (Fig.4(b)).  Our  analytic  derivations  show  that  P=0  at  a=0  is  a  consequence 
of  the  mirror  symmetry  of  the  Taylor  lattice,  which  is  lacking  in  the  dipole  configuration. 
Furthermore,  the  slope  of  the  p(a)  curve  at  o=0  scales  as  n  2,  where  n  is  the  number  of 
dislocations  in  the  Taylor  lattice.  This  means  that  one  cannot  assume  the  dislocation 
multipole  structure  is  a  simple  sum  of  dipoles,  each  making  an  identical  contribution  to  p. 

The  finding  from  this  analysis  completely  changes  our  understanding  of  the  dislocation 
mechanism  to  acoustic  nonlinearity  in  fatigued  metals.  We  expect  them  to  have  profound 
impact  on  the  development  of  the  ultrasonic  non-dcstructive  evaluation  method  for  metal 
fatigue  based  on  the  acoustic  nonlinearity  effect. 


3.  Dislocations  in  Anisotropic  Elastic  Media 

Most  of  the  engineering  metals  are  elastically  anisotropic.  For  simplicity,  most  of  the 
dislocation  dynamics  simulation  programs,  including  ParaDiS  (used  in  this  project), 
assume  isotropic  elasticity.  We  have  designed  and  implemented  efficient  algorithms  to 
enable  dislocation  dynamics  simulations  under  anisotropic  elasticity.  The  algorithms 
have  been  implemented  and  extensively  tested  in  a  Matlab  code  (DDLab)  and  has 
recently  been  implemented  in  ParaDiS  [6].  The  resulting  algorithm  is  significantly  faster 
than  previous  attempts  (e.g.  Rhee  et  al.  Mater.  Sci.  Eng.  A  309-310,  288,  2001). 

4.  Collaboration  with  AFRL  scientists: 

I  am  in  contact  with  Satish  Rao,  Mike  Uchic,  Dennis  Dimiduk,  Chris  Woodward  in 
Wright-Patterson  AFB,  helping  them  with  the  ParaDiS  program.  I  have  visited  these 
scientists  in  WPAFB  in  Feb.  2009  to  discuss  research  collaborations  on  dislocation 
dynamics.  I  have  also  discussed  with  Satish  Rao  and  Mike  Uchic  during  the  annual 
grantee’s  meeting  in  July  2009  in  Fairborn,  OH. 
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